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Abstract 



A measurement of the single-top-quark f-channel production cross section in pp col- 
lisions at \fs = 7 TeV with the CMS detector at the LHC is presented. Two different 
and complementary approaches have been followed. The first approach exploits the 
distributions of the pseudorapidity of the recoil jet and reconstructed top-quark mass 
using background estimates determined from control samples in data. The second 
approach is based on multivariate analysis techniques that probe the compatibility 
of the candidate events with the signal. Data have been collected for the muon and 
electron final states, corresponding to integrated luminosities of 1.17 and 1.56 fb -1 , re- 
spectively. The single-top-quark production cross section in the f-channel is measured 
to be 67.2 ± 6.1pb, in agreement with the approximate next-to-next-to-leading-order 
standard model prediction. Using the standard model assumption that the CKM ma- 
trix element |V t t,| is less than one, we infer the interval 0.92 < |y t t,| < 1 at the 95% 
confidence level. 
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1 Introduction 

Single top quarks can be produced through charged-current electroweak interactions. Due to 
the large top-quark mass, these processes are well suited to test the predictions of the standard 
model (SM) of particle physics and to search for new phenomena. Measurements of the single- 
top-quark production cross section also provide an unbiased determination of the magnitude 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element | V t b I • 

Single-top-quark production was observed in proton-antiproton collisions at the Tevatron col- 
lider with a centre-of-mass energy of 1.96 TeV |Q]-[3]|. The cross section increases by a factor of 
20 at the Large Hadron Collider (LHC) with respect to the Tevatron. The first measurements 
of the single-top-quark production cross section in proton-proton collisions at a centre-of-mass 
energy of 7 TeV were performed by the Compact Muon Solenoid (CMS) [4] and ATLAS HI 
experiments. 

Previous measurements are compatible with expectations based on approximate next-to-leading- 
order and next-to-next-to-leading logarithm (NLO+NNLL) perturbative quantum chromody- 
namics (QCD) calculations. In these calculations, three types of parton scattering processes are 
considered: f-channel and s-channel processes, and W-associated single-top-quark production 
(tW). The dominant contribution to the cross section is expected to be from the f-channel pro- 
cess with a cross section of = 64.6+Qy ~^\y pb for a top-quark mass of m t = 172.5 GeV/c 2 . 




Figure 1: Leading order Feynman diagrams for single-top-quark production in the f-channel: 
2^2 (left) and 2 -> 3 (right). 

This paper extends the previous CMS measurement [4] of the f-channel cross section. The 
single-top-quark production cross section measurement is based on pp collision data at y/s = 
7 TeV collected during 2011 with the CMS experiment, corresponding to integrated luminosities 
of 1.17 and 1.56 ftT 1 with muon and electron final states, respectively. Events with leptonically 
decaying W bosons are selected: t — > bW — > b£v (£ = e or This measurement is used to 
determine the CKM matrix element | V t b I • 

The f-channel event signature (Fig. [lj typically comprises one forward jet scattered off a top 
quark. The decay products of the top quark mainly appear in the central region of the detec- 
tor. A dedicated event selection is applied, and then measurements with two complementary 
approaches are performed. The first approach exploits the reconstructed top-quark mass and 
one of the angular properties specific to f-channel top-quark production: the forward pseudo- 
rapidity distribution of the light jet recoiling against the top quark. This analysis is referred 
to as the \r]ji \ analysis. It is straightforward and robust and has little model dependence. The 
second approach exploits, via multivariate discriminators, the compatibility of the signal can- 
didates with the event characteristics predicted by the SM for electroweak top-quark produc- 
tion. This approach aims for a precise f-channel cross section measurement by optimising the 
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discrimination between signal and background. The systematic uncertainties are constrained 
by simultaneously analysing phase space regions with substantial i-channel single-top-quark 
contributions and regions where they are negligible. Because of the higher complexity of this 
approach, two independent multivariate analyses are conducted that cross-check each other. 
One is based on a Neural Network (NN) and the other on Boosted Decision Trees (BDT). After 
validating the consistency of the results, the final result is determined by combining the three 
analyses using the Best Linear Unbiased Estimator (BLUE) QQ. 

2 CMS Detector 

The CMS Detector is described in detail elsewhere H), an d the key components for this anal- 
ysis are listed below. The main feature of the CMS detector is a superconducting solenoidal 
magnet with a length of 13 m and a diameter of 6 m which provides an axial magnetic field of 
3.8 T. The bore of the solenoid contains various particle detection systems. Charged particle 
trajectories are recorded by the silicon pixel and strip tracker covering the region \t]\ < 2.5, 
where the pseudorapidity rj is defined as t] = — ln[tan(0/2)], with 9 the polar angle of the 
trajectory of the particle with respect to the anticlockwise beam direction. A crystal electro- 
magnetic calorimeter (ECAL) and a brass/ scintillator hadronic calorimeter (HCAL) surround 
the tracking volume and extend to \t] \ < 3.0. A quartz-fibre and steel absorber Cherenkov-light 
hadron forward detector extends the calorimeter coverage to \tj\ < 5.0. Muons are detected in 
gas-ionization detectors embedded in the steel return yoke outside the solenoid. The detector 
is nearly hermetic, which permits good energy balance measurements in the plane transverse 
to the beam line. 

3 Collision Data and Simulation 

The present analyses use data sets collected during 2011 corresponding to integrated luminosi- 
ties of 1.17 and 1.56 fb _1 for the muon and electron channels, respectively. 

Single-top-quark f -channel signal events from Monte Carlo (MC) simulation used in this study 
have been generated with the next-to-leading-order (NLO) MC generator POWHEG |flQHT3l in- 
terfaced to PYTHIA 6.4.24 [14] for the parton showering simulation. This ensures good mod- 
elling of the 2^2 (Fig. [1} left) process and all NLO contributions, in particular of the 2— >3 
diagram (Fig.JTJ right). Signal events have also been generated using COMPHEP [15| to study 
systematic uncertainties related to the MC generator. A top-quark mass of 172.5 GeV/c 2 has 
been assumed in the simulations. 

Several SM processes are taken into account as background to the analysis. The POWHEG gen- 
erator interfaced with PYTHIA is used to model the single-top-quark tW and s-channel events, 
which are considered as background in this analysis. The tree-level matrix-element generator 
MADGRAPH [16] interfaced to PYTHIA is used for top-quark pair (tf) production and for inclu- 
sive single boson production V + X, where V = W or Z, and X can indicate one or more light or 
heavy partons. In the following, we will explicitly distinguish the processes with V+light par- 
tons, Wc, and VQQ, with Q = b or c. A procedure implemented during the event generation 
and based on the "MLM prescription" fFT) avoids double counting between matrix element 
and parton shower generated jets. The remaining background samples are simulated using 
PYTHIA. These include diboson production (WW, WZ, ZZ), 7+jets, and multijet QCD enriched 
events with electrons or muons coming from the decays of b and c quarks, as well as muons 
from the decay of long-lived hadrons. The CTEQ6 |fl8| parton distribution functions are used 
for all simulated samples. 
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All generated events undergo a full simulation of the detector response according to the CMS 
implementation of Geant4 |19| , and are processed by the same reconstruction software used 
for collision data. 

4 Event Selection and Reconstruction 

Events are characterized by a single isolated muon or electron and momentum imbalance due 
to the presence of a neutrino, with one central b jet from the top-quark decay. An additional 
light-quark jet from the hard-scattering process is often present in the forward direction. A 
second b jet produced in association with the top quark can also be present (Fig. [IJ right), 
although it yields a softer pj spectrum with respect to the b jet coming from the top-quark 
decay. 

The trigger used for the online selection of the analysed data for the muon channel is based on 
the presence of at least one isolated muon with a transverse momentum pj > 17GeV/c. For the 
electron channel, an isolated electron trigger with a transverse momentum pj > 27GeV/c was 
used for the initial data-taking period, corresponding to an integrated luminosity of 216 pb _1 . 
For the remaining data-taking period, a trigger selecting at least one electron with pj > 25 GeV/ c 
and a jet with pj > 30 GeV/c was used. The jet is identified in the trigger processing as coming 
from a fragmentation of a b quark using the Track Counting High-Efficiency (TCHE) b-tagging 
algorithm described in Ref. |20| . 

At least one primary vertex is required, reconstructed from a minimum of four tracks, with a 
longitudinal distance of less than 24 cm, and a radial distance of less than 2 cm from the inter- 
action point. To select good muon and electron candidates, the same lepton identification as 
described in Ref. [21] is applied. Electrons, muons, photons, and charged and neutral hadron 
candidates are reconstructed and identified using the CMS particle-flow (PF) algorithm [22 J. 
The missing transverse momentum vector tfj is reconstructed from the momentum imbalance 
of PF particle candidates in the plane transverse to the beams direction. The magnitude of 
is the missing transverse energy tfj. The presence of exactly one isolated muon or electron 
candidate originating from the primary vertex is required in the event. Muon candidates are 
selected by requiring a transverse momentum pj > 20GeV/c and a pseudorapidity \rj\ < 2.1. 
Electron candidates must have a transverse momentum pj > 30GeV/c and \rj\ < 2.5. Lep- 
ton isolation is defined as the sum of the transverse energy deposited by stable charged 
hadrons, neutral hadrons, and photons in a cone of AR = a/(A?/) 2 + (A<p) 2 = 0.4 around the 
charged-lepton track, divided by the transverse momentum of the lepton. Muon isolation is 
ensured by requiring 1^ < 0.15, while the isolation requirement is tightened to 0.125 for elec- 
trons. An electron candidate is rejected if it is identified as originating from the conversion 
of a photon into an electron-positron pair or if it fails the identification criteria described in 
Ref. ||2T]. Events are also rejected if an additional muon candidate is present that passes looser 
quality criteria, namely, pj > 10 GeV/c, \rj\ < 2.5, and I^ el < 0.2. For additional electrons, the 
required transverse momentum is pj > 15GeV/c. 

Jets are defined by clustering PF candidates according to the anti-fcx algorithm fl23"| with a dis- 
tance parameter of 0.5. The analysis considers jets whose calibrated transverse momentum is 
greater than 30GeV/c for \rj\ < 4.5. An event is accepted for further analysis only if at least 
two such jets are reconstructed. A jet is identified as coming from a b-quark fragmentation if 
it passes a tight threshold on the Track Counting High-Purity (TCHP) b-tagging algorithm 1120 1 
corresponding to a misidentification probability of 0.1%. The difference between simulated 
and measured b-tagging efficiencies, for true and misidentified b jets, is corrected by scaling 
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the simulated event yields according to px-dependent scale factors determined from control 
samples from the data fl20|. 

Events with a muon that is not from a leptonic decay of a W boson are suppressed by requiring 
a reconstructed transverse W boson mass ?«x = W2p T # T (l - cos(A<^£ )) > 40GeV/c 2 , where 

A(p e is the azimuthal angle between the muon and the tfj directions and px is the transverse 
momentum of the muon. For the electron channel, where the QCD multijet contamination is 
larger, the requirement %i > 35 GeV is applied instead of the mj selection. 

To classify signal and control samples, different event categories are defined and denoted as 
"n-jets m-btags", where n is the number of selected jets (2, 3, or 4) and m is the number of se- 
lected b-tagged jets (0, 1, or > 2). The single-top-quark f-channel signal is primarily contained 
in the category "2-jets 1-btag", followed by "3-jets 1-btag", as the second b jet, which is pro- 
duced in association with the top quark, is mostly out of acceptance. The other categories are 
dominated by background processes with different compositions. In particular, the "2/3-jets 
0-btags" categories are enriched in events with a W boson produced in association with light 
partons (u, d, s, g). The "3-jets 2-btags" and "4-jets 0/1/2-btags" categories are enriched in tt 
events. 

To extract the signal content, the NN and BDT analyses utilize the following six categories si- 
multaneously for the measurement of the signal cross section: "2-jets 1-btag", "3-jets 1-btag", 
"4-jets 1-btag", "2-jets 2-btags", "3-jets 2-btags", and "4-jets 2-btags". The latter four categories 
are used to constrain nuisance parameters (e.g. the b-tagging efficiency or background normal- 
isation). Besides these six categories, the three 0-tag categories are used to check the modelling 
of input variables. The (77,7 1 analysis extracts the signal content from the "2-jets 1-btag" cate- 
gory only, but uses the "2-jets 0-btags" and "3-jets 2-btags" categories to check the modelling 
of backgrounds. 

The reconstruction of the top quark from its decay products leads to multiple choices of pos- 
sible top-quark candidates. In the first step, the W-boson candidate is reconstructed from the 
charged lepton and from tfj following the procedure described in Ref. 0|. In the second step, 
the top-quark candidate is reconstructed by combining the W-boson candidate with a jet identi- 
fied as coming from a b quark. Depending on the analysis category, the ambiguity in the choice 
of the b-quark jet from the top-quark decay and the recoiling light quark has to be resolved. 
Events in the "2-jets 1-btag" category have no ambiguity: the b-tagged jet is associated with the 
top-quark decay, and the other jet is considered to be a light-quark jet. 

In the other categories, the top quark and light-quark jet reconstruction has been optimized 
for the purpose of each analysis and differs among them. In the NN and BDT analyses, the 
most forward jet is chosen to be the light-quark jet in categories where two jets and zero or 
two b tags are required. The other jet is associated with the top-quark decay. In categories 
where three or more jets are required, in the case of one or two required b tags, the jet with the 
lowest value of the TCHP discriminator is assumed to originate from the light quark. If no b 
tag is required, the most forward jet is associated with the light quark. From the remaining jets, 
the one which together with the reconstructed W boson has a reconstructed top-quark mass 
closest to m t = 172.5 GeV/c 2 is chosen as the b jet coming from the top-quark decay. In 
the I rjji I analysis, the jet with the highest value of the TCHP discriminator is used for the top- 
quark reconstruction in the "2-jets 0-btags" and "3-jets 2-btags" categories. The inclusive \rj\ 
distribution of both jets is used in "2-jets 0-btags", while for the "3-jets 2-btags" category the 
I ?/ 1 of the non-b-tagged jet is used. 
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In the | rjji | analysis the invariant mass mg v \, of the reconstructed top quark is used to further 
divide the "2-jets 1-btag" category into a f-channel enriched signal region (SR), defined by 
selecting events within the mass range 130 < m ( v \, < 220 GeV/ c 2 , and a W boson and tt enriched 
sideband region (SB), defined by selecting events that are outside this m ^ mass window. The 
event yield in the SR is summarised in Table [T] for the muon and electron channels, together 
with expectations from simulated signal and backgrounds, and for QCD multijet events, which 
are determined from control samples of data, as described in Section [^Tj 



Table 1: Event yield with statistical uncertainties of the \rjj'\ analysis for the signal and main 
background processes in the signal region, after applying the mass requirement for the \i 
and e channels. The yields are taken from simulation except for the QCD multijet yield, which 



is obtained from control samples of data as described in Section 5.1 The normalisation of the 



W + jets processes obtained from control samples of data is further discussed in Section 5.2 



Process 


Muon yield 


Electron yield 


f-channel 


617 ±3 


337 ± 2 


tW channel 


107 ± 1 


70.2 ± 0.9 


s-channel 


25.6 ± 0.5 


14.7 ± 0.4 


tt 


661 ±6 


484 ±5 


W + light partons 


92 ±7 


38 ±4 


Wc(c) 


432 ± 14 


201 ± 9 


Wb(b) 


504 ± 14 


236 ± 10 


Z + jets 


87 ±3 


13 ± 1 


Dibosons 


23.3 ± 0.4 


10.7 ± 0.3 


QCD multijet 


77 ±3 


62 ±3 


Total 


2626 ± 22 


1468 ± 16 


Data 


3076 


1588 



5 Background Estimation and Control Samples 

Several control data samples are used for two main purposes: 

• to check the distributions of variables used as inputs to the analyses and the agree- 
ment between data and simulation; 

• to determine from data the yields and distributions of variables of interest for the 
main background processes. 

All the analyses use rate and shape determinations from QCD multijet background data. The 
| rjji | analysis also determines the yields and distributions of the background processes by us- 
ing W boson production in association with jets from light quarks as well as c and b quarks. 
The NN and BDT analyses take the shape of the W+jets background from MC simulation, but 
consider the impact of shape uncertainties in the evaluation of the systematic uncertainties. 
In all three analyses the rate of the W+jets background is determined "in situ" by the signal 
extraction method as described in Section |U 

5.1 QCD Multijet Background Estimation 

The yield of the QCD multijet background in the different categories is measured by perform- 
ing fits to the distributions of the transverse W boson mass nij in the muon channel and Ut in 
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the electron channel. A maximum-likelihood fit to the distribution of mj or 5t is performed 
assuming the following parameterisation: F(x) = a ■ S(x) + b ■ B(x), where x is mj or for 
the muon and electron channels, respectively, and S(x) and B(x) are the expected distributions 
for the sum of all processes including a W boson in the final state and for QCD multijet events, 
respectively. Function S(x) is taken from simulation, while B(x) is extracted directly from data. 
The parameters a and b are determined from the fit. The QCD multijet background yield is es- 
timated to be the area of the fitted curve b ■ B(x) in the range mj > 40GeV/c 2 for the muon 
channel and %j > 35 GeV for the electron channel, as mentioned in Section |4| 

QCD multijet enriched samples from data are used to model the distributions B(x) of mj and 
For the muon channel, this sample is selected by inverting the isolation requirement. For 
the electron channel, the selected electron is required to fail at least two of the three following 
quality requirements: < 0.1, the distance of closest approach to the primary vertex on the 
x-y plane 8 X y < 0.02 cm, and the electron identification criteria given in Ref. (2TJ. It was verified 
by simulation that the mj and distributions for QCD multijet like events are not significantly 
affected by this altered event selection. 

In the | t]j' | analysis the fits to the mj and distributions cannot be performed reliably in the 
SR and SB separately due to the limited size of the simulated samples, which would introduce 
large uncertainties in the signal modelling. For this reason, the fit is performed on the entire 
"2-jets 1-btag" sample. The number of QCD multijet events in the SB and SR regions is deter- 
mined by scaling the total QCD multijet yield, obtained from the fit, by the fraction of events 
in the two regions (SB and SR) of the m^b distribution, as determined from the QCD multijet 
enriched region. In the simulation, the distributions of the relevant variables obtained from the 
QCD multijet enriched sample are consistent with the ones in the SR. The QCD multijet yields 
restricted to the SR are reported in Table [l] 

For all three analyses, the relative uncertainties on the QCD multijet yield estimates are taken to 
be ±50% for the muon channel and ±100% for the electron channel. Several cross checks have 
been performed; for example, the same fits have been repeated taking B(x) from simulation 
for both channels. Moreover, the choice of x = mj or Uj has been inverted for the muon and 
electron channels, in this way performing fits on mj for the electron channel and for the 
muon channel. The results in each case are in agreement with the previous estimate within the 
assumed uncertainties. 

5.2 W+jets Background Estimation and Other Control Samples 

A check of the modelling for the W+jets background is carried out for each of the three analyes 
in the 0-tag control regions. In particular, the "2-jets 0-btags" category is highly enriched in 
W+ light jet events. The modelling of tt background is checked in the "3-jets 2-btags" as well as 
the "4-jets" categories. In general, the event yields are within the systematic uncertainties are 
reasonably well reproduced by the simulation. The shapes of the relevant variables, | t]f | and 
ffiMv an d the input variables of the NN and BDT analyses show good agreement between data 
and simulation. 

Table [Ilshows a difference between the total observed and expected yields for the \t]f | analysis. 
This difference can be attributed to excesses in data for the Wb+X and Wc+X processes. The 
ATLAS collaboration reported in Ref. [l24l| that the fiducial W+b-jet cross section in the lepton 
and one or two jets final state is a factor of 2.1 larger than the NLO prediction, but is still 
consistent at the level of 1.5 standard deviations with this SM prediction. 

Motivated by these observed excesses in comparison to the SM NLO calculations, the | rjji | anal- 



7 



ysis determines the W+jets background yield and | rjji | distribution from data. The | rjy | distribu- 
tion for W+jets process is extracted from the SB by subtracting the | r\v | distribution of all other 
processes from the data. The event yield and \r\v \ distributions used for these subtractions are 
taken from simulations of tt, single-top-quark s- and tW- channels, and diboson production. 
The QCD multijet event yield an d \ru > | distribution are extracted from data and extrapolated 
to the SB as described in Section 5.1 The \rjA distribution for W+jets processes in the SB is 
therefore used in the SR for the signal extraction procedure (see Section [6]), assuming that the 
shapes in the SB and SR are compatible with each other. For the muon channel, the compati- 
bility of the distributions in the two regions has been verified through a Kolmogorov-Smirnov 
compatibility test, yielding a p-value of 0.47, and a % 2 test, yielding a p-value of 0.63. For the 
electron channel, the Kolmogorov-Smirnov compatibility test has a p-value of 0.51 and the x 2 
test a p-value of 0.60. The stability of the extracted shape has been tested by varying the sample 
composition in terms of tt and signal fractions by 20% and 100%, respectively. The extracted 
shapes are compatible with a p-value greater than 0.9 in both cases. 



6 The \tjj/\ Analysis 

The signal yield is extracted using a maximum-likelihood fit to the observed distribution of 
| rjji | . The signal distribution for the fit is taken from simulation. The W/ Z+jets component of 
the background is normalised to the value obtained from the extraction procedure described in 



Section 5.2 and then added to the diboson processes, resulting in the electroweak component 
of the background for the fit. The signal and the electroweak components are unconstrained 



in the fit, whereas the QCD multijet component is fixed to the result determined in Section 5.1 
A Gaussian constraint is applied to tt and other top-quark backgrounds. Figure [2] shows the 
distribution of \r]ji \ obtained from the fit. Figure [3] shows the distribution of the reconstructed 
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Figure 2: Result of a simultaneous fit to \rjA in the muon (left) and electron (right) decay 
channels. 



top-quark mass ni( v \, normalised to the fit results, restricting to the highly enriched region of 
single-top-quark events for \tjji\ > 2.8. 

Single-top-quark i-channel production at the LHC is expected to be characterised by two fea- 
tures. First, the top-quark cross section is about a factor two larger than the top-anti-quark cross 
section 0. This can be experimentally accessed via the charge of the muon or electron. Second, 
top quarks are almost 100% polarised with respect to a certain spin axis due to the V-A nature 
of the couplings. This can be studied via the cos 8* distribution II25II , where 9* is defined as the 
angle between the charged lepton and the non-b-tagged jet, in the reconstructed top-quark rest 
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Figure 3: Distributions of requiring \r\ii\ > 2.8, for muons (left) and electrons (right), ob- 
tained by normalising each process yield to the value from the fit. Because of limited simulated 
data, the background distribution is smoothed by using a simple spline curve. 
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Figure 4: Distinct single-top-quark i-channel features in the SR for \rjv \ > 2.8, for the electron 
and muon final states combined. The charge of the lepton (left) and cos 8* (right). All processes 
are normalised to the fit results. Because of limited simulated data, the background distribution 
is smoothed by using a simple spline curve (right). 

frame. The observed charge asymmetry and the cos0* distribution are presented in Fig. |4] for 
muon+electron events in the SR, for \rjj> \ > 2.8. 

7 Neural Network Analysis 

In the NN analysis, several kinematic variables, which are characteristic of SM single-top-quark 
production, are combined into a single discriminant by applying an NN technique. The NEU- 
ROBAYES package |26l|2"7] used for this NN analysis combines a three-layer feed-forward NN 
with a complex, but robust, preprocessing. To reduce the influence of long tails in distributions, 
input variables are transformed to be Gaussian distributed. In addition, a diagonalisation and 
rotation transformation is performed such that the covariance matrix of the transformed vari- 
ables becomes a unit matrix. To obtain good performance and to avoid overtraining, the NN 
uses Bayesian regularisation techniques for the training process. The network input layer con- 
sists of one input node for each input variable plus one bias node. The hidden layer is adapted 
to this particular analysis and consists of one more node than the input layer. The output node 
gives a continuous discriminator output in the interval [—1, 1]. For the training of the NN, 
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after applying the full event selection as described in Section |4| simulated samples of signal 
i-channel single-top-quark events and background events from tt, W+jets, and Z+jets samples 
are used. The ratio of signal to background events in the training is chosen to be 50:50 and the 
background processes are weighted according to the SM prediction, as outlined in Section [5] 
The NN is trained such that i-channel single-top-quark events tend to have discriminator val- 
ues close to 1, while background events tend to have discriminator values near —1. Because 
of their different event selections, separate neural networks are trained for muon and electron 
events. 

During the preprocessing, the training variables are ranked by the significance of their correla- 
tion to the target discriminator output. The correlation matrix of all training variables and the 
target value is calculated. The variable with the smallest correlation to the target is removed 
and the loss of correlation is calculated. This is repeated until the correlation of all variables 
with the target is determined. The significance of a variable is calculated by dividing the loss 
of correlation with the target by the square root of the sample size. In order to select variables 
which contain information that is not already incorporated by other variables, a selection cri- 
terion of > 3 a on the significance has been chosen. A set of 37 variables remains for the muon 
channel when applying this selection criterion. For the electron channel, a set of 38 variables 
remains. The validity of the description of these input variables and the output of the NN dis- 
criminant is confirmed in data with negligible signal contribution. Furthermore, it is verified 
with a bootstrapping technique [28] that the bias of the cross section measurement due to a 
possible overtraining of the NN is negligible. 



CMS ,s = 7TeV L= 1.17 lb' 



Data 

1-channel 

It, s-channel, tW 

W/Z + jets, Diboson 

QCD multijet 




5 700 
> 
LU 

600 



500 
400 
300 



200 
100 



CMS ,s = 7TeV L= 1.17 tb 



Muon, "3-jets 1-btag" 



Data 
^| t-channel 

tt, s-channel, tW 
^| W/Z + jets, Diboson 

QCD multijet 



-1 -0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 1 
NN output 



-1 -0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 1 
NN output 



Figure 5: Distributions of the NN discriminator output in the muon channel for the "2-jets 
1-btag" (left) and "3-jets 1-btag" (right) categories. Simulated signal and background contribu- 
tions are scaled to the best fit results. 

The variables with the highest ranking in both networks are \rjj> |, mj, the invariant mass of the 
two leading jets, and the total transverse energy of the event. 

The distributions of the NN discriminator in the "2-jets 1-btag" and "3-jets 1-btag" signal cate- 
gories are shown in Fig. |5]for the muon channel and in Fig. |6]for the electron channel. For the 
remaining categories "2-jets 2-btag", "3-jets 2-btags", "4-jets 1-btag", and "4-jets 2-btags", the 
events are combined into one plot in Fig.[7|for the muon channel on the left and for the electron 
channel on the right. 



8 Boosted Decision Trees Method 



The BDT method was previously used in the first CMS measurement of the i-channel cross 
section (4J. The current analysis uses a significantly larger data sample. It has been optimised 
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Figure 6: Distributions of the NN discriminator output in the electron channel for the "2-jets 
1-btag" (left) and "3-jets 1-btag" (right) categories. Simulated signal and background contribu- 
tions are scaled to the best fit results. 
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Figure 7: Distributions of the NN discriminator output in the background dominated region. 
All events from the signal depleted categories "2-jets 2-btags", "3-jets 2-btags", "4-jets 1-btag", 
and "4-jets 2-btags" are combined for the muon channel (left) and the electron channel (right). 
Simulated signal and background contributions are scaled to the best fit results. 
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using a "blind" analysis with the optimisation based exclusively on regions where the signal 
contribution is negligible and all the selections frozen before the signal region has been looked 
at. It further provides an increase of the measurement sensitivity and a reduction of systematic 
uncertainties. The BDT analysis was designed following Ref. |29|. 

The adopted BDT algorithm constructs 400 decision trees using the Adaptive Boosting algo- 
rithm as implemented in Ref. ||30l . The BDT training is carried out separately for the electron 
and muon final states, for each of the single-top-quark "2-jets 1-btag" and "3-jets 1-btag" signal- 
enriched regions, to provide a total of four BDTs. The background processes are input to the 
training according to their theoretical cross section predictions as outlined in Section [5] One 
third of the simulated signal and background samples are used for the training, one third are 
used to verify the performance of the trained BDT, while the remaining simulated data provide 
an unbiased sample used for the cross section evaluation. 

Eleven observables reconstructed in the detector are chosen based on their power to discrimi- 
nate between signal and background events. 

The adopted variables are the lepton transverse momentum; the pseudorapidities of the most 
forward non-b-tagged jet and of the jet with the highest transverse momentum; the invariant 
mass of all reconstructed jets in the event; the angular separation between the two leading 
jets; the sums of the hadronic energy and of the hadronic transverse energy; the reconstructed 
top-quark mass using the jet with the highest b-tag discriminator; the reconstructed top-quark 
mass using the jet giving the reconstructed top-quark mass closest to 172GeV/c 2 ; the cosine 
of the angle between the reconstructed W boson, in the rest frame of the sum of four-vectors 
of the W boson and leading jet, and the sum of four-vectors of the W boson and leading jet; 
and the sphericity of the event. The validity of the description of these input variables and the 
output of the BDT classifiers is confirmed in data with negligible signal contribution using a 
Kolmogorov-Smirnov test. 

The QCD multijet background evaluation is described in Section |5j The determination of the 
single-top-quark production cross section in the i-channel, including the treatment of statistical 
and systematic uncertainties, is performed by using the classifier distributions in all 12 analysis 
categories simultaneously. 

The measured distributions of the classifier outputs in the "2-jets 1-btag" and "3-jets 1-btag" 
signal enriched categories are shown separately for the muon and electron channel in Figs. [8] 
and [9] The measured distributions of the classifier outputs in the "2-jets 2-btags","3-jets 2- 
btags", "4-jets 1-btag", and "4-jets 2-btags" signal depleted categories are shown separately for 



the muon and electron channels in Fig. 10 



9 Determination of the Cross Section with Multivariate Analyses 

The signal cross section is determined in the NN and BDT analyses from the data distribu- 
tions of the corresponding multivariate discriminator, modelled as 1-dimensional histograms, 
in the different categories, which are used simultaneously. The distributions for signal and 
backgrounds are taken from simulation. The signal yield is measured in terms of the signal 
strength ]i, defined as the actual cross section divided by the SM prediction. A binned likeli- 
hood function based on a Poisson distribution is then constructed from the data. The simulated 
values are taken as the Poisson mean values. This likelihood function p(data|p) is related to 
the posterior distribution p(p|data) of p using Bayes' theorem, up to a normalisation factor 



p(p|data) <x p(data|p) • 7r(p) 
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Figure 8: Distributions of the BDT discriminator output in the muon channel for the "2-jets 
1-btag" (left) and "3-jets 1-btag" (right) categories. Simulated signal and background contribu- 
tions are scaled to the best fit results. 




BDT output BDT output 

Figure 9: Distributions of the BDT discriminator output in the electron channel for the "2-jets 
1-btag" (left) and "3-jets 1-btag" (right) categories. Simulated signal and background contribu- 
tions are scaled to the best fit results. 




Figure 10: Distributions of the BDT discriminator output in the signal depleted categories for 
the muon channel (left) and electron channel (right). Simulated signal and background contri- 
butions are scaled to the best fit results. 



13 



where rc(]d) denotes a uniform prior distribution of the signal strength. 

To take into account systematic uncertainties, the likelihood definition is initially extended to 
an (m + 1) -dimensional likelihood function which includes the impact of the m relevant sys- 
tematic uncertainty contributions, p'(data|p., 6), where 6 = ■ ■ ,9 m ) contains k normali- 
sation parameters and the m — k nuisance parameters influencing the shape of the simulated 
distributions. The prior distribution of the nuisance parameters n{Q) can be factorised into 
the product of prior distributions Ki{9{) for each nuisance parameter, i = 1, • • • ,k. For nor- 
malisation uncertainties, 7T ; -(0,) are log-normal distributions with average fixed to 1 and given 
standard deviations. For shape nuisance parameters, 7r,(#,) includes the simulated discrimina- 
tor distributions varied by ±1 standard deviation in the corresponding uncertainty where the 
uncertainty is varied according to a Gaussian distribution. To obtain continuous variations in 
each bin, an interpolation algorithm is used to vary histograms between the nominal and the 
±l£T,-varied distributions. 

The (m + 1) -dimensional phase space of the likelihood function takes into account systematic 
uncertainties including their correlations and is constrained by the data. The following poste- 
rior distribution p(f/|data) results, up to a normalisation factor, from integrating the product 
of the likelihood function and the prior distributions in all m dimensions except for the signal 
strength ji: 

pO|data) <x j p'(data|jn,#) • Tt(}i)n(d) d6 . 

The integration is performed using the Markov chain MC method as implemented in the pack- 
age Theta |3H . From this distribution we extract y. as the central value and determine the 
central 68% quantile as the total uncertainty in the measurement including the statistical and 
systematic uncertainties. 

The marginalisation procedure provides an overall uncertainty corresponding to both the finite 
number of events in the analysed sample and the "in situ" constraint of the nuisance param- 
eters. To avoid underestimating the systematic uncertainties, which can arise from the like- 
lihood function used in the Bayesian method by applying too many constraints, checks have 
been performed on the model for experimental uncertainties, and theoretical uncertainties are 
not marginalised. The contribution of those sources of uncertainty that are not marginalised is 
estimated by varying each nuisance parameter by ±1 a separately. The signal extraction proce- 
dure is repeated, with the corresponding change of signal and background shapes. The average 
shift is estimated from many generated pseudo-experiments. An estimate of the statistical un- 
certainty is taken from the total uncertainty after subtracting, in quadrature, the individual 
contributions as evaluated above. 

10 Systematic Uncertainties and Measurement Sensitivity 

For the | rjv | analysis, each systematic uncertainty is evaluated by generating pseudo-experiments, 
which take into account the effect of the corresponding systematic source on the distribution 
of \tjji | and on the event yield of the physics processes. Pseudo-experiments are generated sep- 
arately with templates varied by ±1 a of the corresponding uncertainty. A fit to \t]ji\ is then 
performed on each pseudo-experiment. The mean shift of the fit results, with respect to the 
value obtained in the nominal scenario, is taken as the corresponding uncertainty. 

In the BDT and NN analyses the experimental systematic uncertainties (excluding the lumi- 
nosity) are marginalised with the Bayesian method. Theoretical uncertainties, however, are not 
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10 Systematic Uncertainties and Measurement Sensitivity 



marginalised, but are estimated by generating pseudo-experiments using separate templates, 
varied by ±1 a for the corresponding uncertainty, for each source of systematic uncertainties, 
and repeating the signal extraction procedure. 

A particular uncertainty, which is only present for the \rjv | analysis, concerns the extraction of 
W+jets from data. It is evaluated by generating pseudo-experiments in the SB and repeating the 
signal extraction procedure and fit to | rjj> \ . This method exploits the ansatz that the distribution 
of \t]ji | is the same in both the SR and SB. The uncertainty is taken as the root mean square of 
the distribution of fit results obtained in this way. This uncertainty depends on the amount 
of available data in the SB and is uncorrelated between the muons and electron samples. In 
addition, alternative | tjji | shapes are derived in the simulation by varying the Wb+X and Wc+X 
fractions of the background by factors of ±30% independently in the SR and SB regions. The 
fit procedure is repeated using the new shapes and the maximum difference in the result with 
respect to the central value is added in quadrature to the other uncertainties. 

The following sources of systematic uncertainties are considered in all three analyses. Differ- 
ences in the | t]v | and the two multivariate analyses are remarked upon where relevant: 

• Jet energy scale (JES). All reconstructed jet four-momenta in simulated events are 
simultaneously varied according to the tj and px-dependent uncertainties on the 
jet energy scale [32]. This variation in jet four-momenta is also propagated to #x- 
For multivariate analyses, the complete parameterisation of JES as in [32 [ is consid- 
ered and all parameters are considered as nuisance parameters and included in the 
marginalisation procedure described in Section [9] 

• Jet energy resolution. A smearing is applied to account for the known difference 
in jet energy resolution with respect to data [33 1, increasing or decreasing the extra 
resolution contribution by the uncertainty on the resolution. 

• b tagging. Both b tagging and misidentification efficiencies in the data are estimated 
from control samples [20 J. Scale factors are applied to simulated samples to repro- 
duce the measured efficiencies. The corresponding uncertainties are propagated as 
systematic uncertainties. For multivariate analyses, b tagging average scale factors 
are constrained using the marginalisation procedure described in Section |9j The 
effect of any remaining, unconstrained, b tagging modelling is determined by mod- 
elling possible variations of the b tagging scale factors as a function of pj and \t]\, 
using different degrees (from one to five) of Chebyshev polynomials. The largest 
observed variation with respect to the nominal result is taken as the systematic un- 
certainty, and is found to be negligible. 

• Trigger. Single lepton trigger efficiencies are estimated with a "tag and probe" 
method [34] from Drell-Yan data. The efficiencies of triggers requiring a lepton plus 
a b-tagged jet are parameterised as a function of the jet pj and the value of the TCHP 
b-tag discriminator. The selection efficiencies have been validated using a reference 
trigger. The uncertainties of the parameterisation and an additional flavour depen- 
dency is propagated to the final result. 

• Pileup. The effect of multiple interactions (pileup) is evaluated by reweighting sim- 
ulated samples to reproduce the expected number of pileup interactions in data, 
properly taking into account in-time and out-of-time pileup contributions. The un- 
certainty on the expected number of pileup interactions (5%) is propagated as a sys- 
tematic uncertainty to this measurement. 

• Missing transverse energy. The modelling uncertainty is propagated to the cross 
section measurement. The effect on Ui measurement of unclustered energy deposits 
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in the calorimeters is included. 

• Luminosity. The luminosity is known with a relative uncertainty of ±2.2% [35]. 

• Background normalisation. The uncertainties on the normalisation of each back- 
ground source are listed below. They are propagated as systematic uncertainties 
in the \rjf \ analysis only for dibosons and s- and tW-channel single-top-quark pro- 
cesses. The remaining backgrounds are estimated from data. The uncertainty on tt 
is used as a Gaussian constraint in the signal-extraction fit. In the multivariate anal- 
yses, normalisation uncertainties are accounted for as a prior probability density 
function for the Bayesian inference using a log-normal model. 

• tt: ±15%, based on the statistical uncertainties in Ref. [36 1. 

• Dibosons, single-top-quark s- and tW-channels: ±30%, ±15%, ±13%, re- 
spectively, based on theoretical uncertainties. 

• W/Z+jets: ±100%, ±50%, and ±30% are taken for W+b/c flavour jets, 
W+light flavour jets, and Z+jets, respectively, consistent with previous 
estimates [36J. In the multivariate analyses, the various W+jets processes 
are considered to be uncorrelated, as are the different jet categories, in 
order to avoid too many model assumptions. 

• QCD multijet: the normalisation and the corresponding uncertainty is 
determined from data (see Section p). 

• Limited MC data. The uncertainty due to the limited amount of MC data in the tem- 
plates used for the statistical inferences is determined by using the Barlow-Beeston 
method EZH2Q- 

• Scale uncertainty. The uncertainties on the renormalisation and factorisation scales 
are studied with dedicated single-top-quark and background samples of W+jets, 
Z+jets, and tt events. They are generated by doubling or halving the renormalisation 
and factorisation scale with respect to the nominal value equal to the Q 2 in the hard- 
scattering process. 

• Extra parton modelling (matching). The uncertainty due to extra hard parton ra- 
diation is studied by doubling or halving the threshold for the MLM jet matching 
scheme II39II for W+jets, Z+jets, and tt from its default. 

• Signal generator. The results obtained by using the nominal POWHEG signal sam- 
ples are compared with the result obtained using signal samples generated by 
COMPHEP. In general, the largest model deviations occur in the kinematic distri- 
butions of the spectator b quark [40|. The differences in the transverse momentum 
distribution of the spectator b quark, at the generator level, between 4-flavour and 
5-flavour scheme (FS) POWHEG HP are more than a factor two smaller than the 
differences between POWHEG-5FS and CompHEP over the whole pj range. Pseudo- 
experiments with simulated CompHEP events are generated and the nominal signal 
extraction procedure with POWHEG-5FS templates is repeated. Half of the observed 
shift of the cross section measurement is taken as the systematic uncertainty due to 
signal generator modelling. 

• Parton distribution functions. The uncertainty due to the choice of the parton dis- 
tribution functions (PDF) is estimated using pseudo-experiments, reweighting the 
simulated events with each of the 40 eigenvectors of the CTEQ6 [ 18 1 PDF set and the 
central set of CTEQ10 [42], and repeating the nominal signal extraction procedure. 
For reweighting the simulated events, the lhapdf |[43l[ package is used. 
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10 Systematic Uncertainties and Measurement Sensitivity 



Table 2: Sources of uncertainty on the cross section measurement. 
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Jet energy scale 
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—0.3/ - 


1 QO/ 
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Jet energy resolution 
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Z 

z 

13 


g 


lx/Tiion "Ffi (T <Tt-">t" -U rprn 

IVlUwlL llliiiit;! i 1CLU. 


-2.2/4 


-2.3% 


±1.9% 




Exper 


Electron trigger + reco. 


-0.6/4 


-0.7% 


±1.2% 


-0.8/ ±0.9% 


Hadronic trigger 


-1.3/4 


-1.2% 


±1.5% 


±3.0% 


d> 


Pileup 


-1.0/4 


-0.9% 


±0.4% 


-0.3/ ±0.2% 


rginal 




$t modelling 


-0.0/4 


-0.2% 


±0.2% 


±0.5% 




W+jets 


-2.0/4 


-3.0% 


-3.5/+2.5% 


±5.9% 


(a 
% 


03 
01 


light flavor (u, d, s, g) 


-0.2/4 


-0.3% 


±0.4% 


n/a 




cC 
Sh 


heavy flavor (b, c) 


-1.9/4 


-2.9% 


-3.5/ ±2.5% 


n/a 




bb 


tt 


-0.9/4 


-0.8% 


±1.0% 


±3.3% 




u 


QCD, muon 


±0.{ 


3% 


±1.7% 


±0.9% 




CO 


QCD, electron 


±0.4% 


±0.8% 


-0.4/ ±0.3% 
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±0.3% 


±0.6% 


±0.5% 




Total marginalised uncertainty 
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n/ a 
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-0.0/ ±3.4% 


n/a 


C 

'5b 


<J 
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Scale, t~, s-, tW channels 


-0.4/4 


-1.0% 


±0.2% 


-2.2/ ±2.3% 


t mar 


S-i 

O 
0) 


Matching, tt 


±1.3% 


±0.4% 


±0.4% 


£-channel generator 


±4.2% 


±4.6% 


±2.5% 


o 


H 


PDF 


±1.3% 


±1.3% 


±2.5% 






Total theor. uncertainty 


-6.3/4 


-4.8% 


-4.9/ ±5.9% 


-5.6/ ±4.9% 


Syst. + 1 


leor. + luminosity uncert. 


-8.1/4 


-7.8% 


-8.1/ ±8.4% 


±10.8% 


Total (stat. + syst. + theor. + lum.) 


-10.1/ ±9.5% 


-9.4/ ±10.0% 


±13.8% 



The jet energy scale and jet energy resolution are fully correlated across all samples. The match- 
ing and scale uncertainties are fully correlated between W+jets and Z+jets, but are uncorrelated 
with tt. 

Table [2] summarises the different contributions to the systematic uncertainty on the combined 
(muon and electron) cross section measurement in the three analyses. For the multivariate anal- 
yses, experimental uncertainties and background rates are constrained from the marginalisa- 
tion procedure (except for the luminosity), as described in Section|9] The remaining theoretical 
and luminosity uncertainties are added separately in quadrature to the total uncertainty. 
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11 Results 

11.1 Individual analyses 

The \f]f \ analysis yields the following cross section measurements for the muon and electron 
channels: 

cr t . ch = 73.3 ± 10.4 (stat. + syst. + lum.) ± 4.0 (theor.) pb (muons), 
ff-ch. = 61.6 ± 13.9 (stat. + syst. + lum.) ± 3.5 (theor.) pb (electrons). 

The two measurements are compatible, taking into account correlated and uncorrelated un- 
certainties. The uncorrelated uncertainties include the W+jets and QCD extraction procedure, 
lepton reconstruction and trigger efficiencies, and the hadronic part of the trigger. 

Combining the muon and electron measurements gives 

tff-ch. = 70.0 ± 6.0 (stat.) ± 6.5 (syst.) ± 3.6 (theor.) ± 1.5 (lum.) pb. 

The Bayesian inference is performed with the data samples and the 50% quantile is calculated 
as the best parameter estimate for the signal strength }i. The 84% and 16% quantiles are quoted 
as upper and lower boundaries for the la credible interval. 

The measured single-top-quark i-channel production cross section in the NN analysis is 

Cf-ch. = 69. (stat. + syst. + lum.) ± 3.6 (theor.) pb (muons), 
Cf-ch. = 65.1^g| (stat. + syst. + lum.) ± 3.5 (theor.) pb (electrons). 

The quoted errors are the total uncertainties including marginalised, unmarginalised, and lu- 
minosity uncertainties. The two measurements are compatible within the uncertainties, after 
properly taking into account the correlated contributions. The combination of the muon and 
electron measurements gives 

a t . cK = 68.1 ± 4.1 (stat.) ± 3.4 (syst.)+^ (theor.) ± 1.5 (lum.) pb. 

The measured single-top-quark i-channel production cross section in the BDT analysis is 

^t-ch. = 66.6+™ (stat. + syst. + lum.)!^ (theor.) pb (muons), 
^t-ch. = 66.4+79 (stat. + syst. + lum.)^ (theor.) pb (electrons), 

and their combination gives 

Of-ch. = 66.6 ± 4.0 (stat.) ± 3.3 (syst.)+J| (theor.) ± 1.5 (lum.) pb. 

The results of the three analyses are consistent with the SM prediction. 
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11 Results 



11.2 Combination 

The results of the three analyses are combined using the BLUE method.The statistical cor- 
relation between each pair of measurements is estimated by generating dedicated pseudo- 
experiments. The correlation is 60% between NN and \t]f\, 69% between BDT and \t]f\, and 
74% between NN and BDT. Correlations for the jet energy scale and resolution, b tagging, and 
#x modelling between \r\v | and the two multivariate analyses are expected to be small. This is 
because the determination of the corresponding nuisance parameters, from the marginalisation 
adopted in the BDT and NN analyses, is dominated by in-situ constraints from data samples 
independent of those used to determine uncertainties in the | rjji | analysis. The assumed corre- 
lation for those uncertainties is taken to be 20%. The correlation has, nevertheless, been varied 
from 0% to 50%, with a corresponding variation of the central value by —0.03 pb, and no appre- 
ciable variation has been observed for the combined uncertainty. For trigger uncertainties, the 
correlation between | mi | and the two multivariate analyses is more difficult to ascertain. Vary- 
ing the correlations in the combination from 0% to 100% results in a variation of the central 
value of 0.03 pb, with no appreciable variation of the combined uncertainty. All other uncer- 
tainties are determined mostly from the same data samples used by the two analyses, hence 
100% correlation is assumed. 

The BLUE method is applied iteratively as previously carried out in Ref. |4J. In each iteration, 
the absolute uncertainty is calculated by scaling the relative uncertainties given in Table [2] with 
the combined value from the previous iteration. This is repeated until the combined value 
remains constant. There are no appreciable changes with respect to the non-iterative BLUE 
method. The 0.03 pb variation in the central value, due to changes in correlation coefficients, 
is added in quadrature to the total uncertainty. However, this results in a negligible additional 
contribution. 

The x 2 obtained by the BLUE combination of the three analyses is 0.19, corresponding to a 
p-value of 0.90. The results of the individual analyses are consistent with each other. 

The combined result of the measured single-top-quark f-channel production cross section at 
v /s = 7TeVis 

a t . cK = 67.2 ± 6.1 pb = 67.2 ± 3.7 (stat.) ± 3.0 (syst.) ± 3.5 (theor.) ± 1.5 (lum.) pb. 

11.3 |y t bl Extraction 

The absolute value of the CKM element |V t b| is determined in a similar fashion to Ref. [4], as- 
suming that |y t a| and |V t s| are much smaller than \V tb \, as |V t b| = J cr t . cb / o~^ ch , where cr*^ 
is the SM prediction calculated assuming |V t b| = 1 [TJ. If we take into account the possible 
presence of an anomalous Wtb coupling, this relation modifies taking into account an anoma- 
lous form factor fi v lj44Tl46l , which is not necessarily equal to 1 in beyond-the-SM models. We 
determine 

l/Ly^tbl =J^ = 1.020 ±0.046 (exp.)± 0.017 (theor.), 

V ^f-ch. 

where the first contribution to the uncertainty includes all uncertainties of the cross section 
measurement, and the second the uncertainty on the SM theoretical prediction. From this re- 
sult, the confidence interval for |V t b|/ assuming the constraint |V t b| < 1 and fi v = 1, is deter- 
mined using the unified approach of Feldman and Cousins 147) to be 

0.92 < | V tb | < 1, at the 95% confidence level. 
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12 Conclusions 

The cross section of i-channel single-top-quark production has been measured in pp collisions 
using 2011 data in semileptonic top-quark decay modes with improved precision compared 
to the previous CMS measurement. Two approaches have been adopted. One approach has 
been based on a fit of the characteristic pseudorapidity distribution of the light quark recoiling 
against the single top quark in the f-channel with background determination from data. The 
other has been based on two multivariate discriminators, a Neural Network and Boosted De- 
cision Trees. The multivariate analyses reduce the impact of systematic uncertainties by simul- 
taneously analysing phase space regions with substantial i-channel single-top-quark contribu- 
tions, and regions where they are negligible. The results are all consistent within uncertainties. 
As a consequence, all three analyses have been combined with the Best Linear Unbiased Esti- 
mator method to obtain the final result. 

The combined measurement of the single-top-quark i-channel cross section is 67.2 ± 6.1 pb. 
This is the first measurement with a relative uncertainty below 10%. It is in agreement with the 
approximate NNLO standard model prediction of 64.6^Q7 -\' 7 pb EI- Figure lT] compares this 
measurement with dedicated i-channel cross section measurements at the Tevatron 11481 149| , 
ATLAS 13, and with the QCD expectations computed at NLO with MCFM in the 5-flavour 
scheme f50| and at NLO+NNLL 0. The absolute value of the CKM matrix element V t \, is mea- 
sured to be \ft v V th \ = Y // ^- ch ./< h ch . = 1.020 ± 0.046 (exp.) ± 0.017 (theor.). Assuming f Lv = 1 
and I V t b| < 1, we measure the 95% confidence level interval 0.92 < |y t b| < 1. 




Figure 11: The single-top-quark cross section in the f -channel vs. centre-of-mass energy. The 
error band (width of the curve) of the SM calculation is obtained by varying the top-quark 
mass within its current uncertainty J3|, estimating the PDF uncertainty according to HEPDATA 
recommendations |5Tf , and varying the factorisation and renormalisation scales coherently by 
a factor of two up and down. 
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